Influence of Electron-Phonon Interaction on 
Spin Fluctuation Induced Superconductivity 
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We investigate the interplay of the electron-phonon and the spin fluctuation interaction in the 
superconducting state of YBa2Cu307. The spin fluctuations are described within the nearly anti- 
ferromagnetic Fermi liquid theory and the phonons are treated using a shell model calculation of all 
phonon branches. The electron-phonon coupling is calculated using rigidly displaced ionic potentials 
screened by a background dielectric constant £00 and by holes within the Cu02 planes. We get a 
superconducting state with d a .2_ !/ 2-symmetry whose origin is antiferromagnetic spin fluctuations. 
The phononic contribution of all phonon modes of the system to the d-wave pairing interaction is 
attractive. This gives a positive isotope effect. The size of the isotope exponent depends strongly 
on the relative strength of the electron-phonon and spin fluctuation coupling. Due to the strong 
electronic correlations no phononic induced superconducting state of s-wave character, is possible. 



I. INTRODUCTION 



The mechanism responsible for the superconducting 
state of high-T c superconductors is still discussed highly 
intensiwsly. Cooper-pairing induced by spin fluctua- 
tionscrEl seems to n be likely. Based on the analysis of 
NMR experiments^, it offers an explanation for transport 
and optical experimc»teE3 and a variety of anomalies in 
underdoped cupratesoQ and furthermore yields a d x 2_ y 2- 
wave pairing state with a transition temperature T c of 
correct order of magnitudes 

However, there are also observations, which indicate ef- 
fects due to electron-phonon interaction. Raman scatter- 
ing experiments and IR measurements show the softening 
of some ghonon modes as the temperature is decreased 
below TJtEj. In YBa2Cu307 the largest shift was found 
for the 340cm -1 (42meV) Raman-active mode of Big- 
symmetry, which involves the out-of-plane vibration of 
the oxaigen atoms 0(H) and O(III) in the copper-oxygen 
planesEx For some oxygen-related phonon modes the 
onset of the frequency softening was found to start even 
50-60 K above T C E3. The temperature dependence of this 
anomaly was found to be similar to that of the NMR re- 
laxation rate (T1T) -1 , which, in the latter case, is due 
to an anomalous behavior of the spin excitations, and 
indicates a possibk-coupling between phonons and spin 
degrees of freedomtd Especially interesting is the obser- 
vation of the anomalous oxygen isotope effect in the high- 
Tc superconductorstij. At optimal doping with largest T c 
only a very small oxygen isotope exponent a = — 0.07 
is observed, whereas in the underdoped and overdoped 
systems the isotope exponent can reach values close to 
the BCS- value a = \. 

Due to these experimental findings it is of interest to in- 
vestigate in more detail the influence of electron-phonon 
interaction on the spin fluctuation induced superconduct- 
ing state. In this context it was shown that oxygen buck- 
ling modes support a d K 2_„2-pairing state, whereas oxy- 



gen breathing modes will support an s-wave stateBH. 
Hence, an analysis of the role of the electron-phonon in- 
teraction in the cuprates has to take into account all vi- 
brational modes. 

In this paper we investigate the interplay of the spin 
fluctuation and the electron-phonon interaction in the su- 
perconducting state of YBa2Cu30y. We find that taking 
into account both spin fluctuation and electron-phonon 
interaction gives a superconducting state of d x 2_ y 2- 
symmetry, while electron-phonon interaction alone would 
yield an anisotropic s-wave pairing with a much smaller 
critical temperature. The d^.^-superconducting state 
is caused by antiferromagnetic spin fluctuations and the 
critical temperature T c is reduced by simoultaneous con- 
sideration of phonons. Nevertheless we find a positive 
isotope effect. We show that the momentum depen- 
dence of the total electron-phonon interaction amplifies 
the pairing interaction by an amount of about 5-10% 
in the d x 2_ y 2--w&ve channel, which is a necessary con- 
dition for a positive isotope effect. Furthermore we iden- 
tify the phonon-modes, which make the most dominant 
contribution to the d^-j^ -pairing interaction. The de- 
crease of T c is explained by the large isotropic part of the 
electron-phonon interaction, which enhances the electron 
scattering rate significantly. However, this effect will be 
reduced when correlation effects are taken into account, 
since correlation induced suppression of charge fluctua- 
tions will predominantly reduce the isotropic part of the 
electron-phonon interaction. Our results are obtained 
with a self-consistent strong coupling Eliashberg calcu- 
lation. Finally, in order to make this interplay of spin 
fluctuation and electron-phonon interaction more trans- 
parent, we present a simple weak coupling model. 

In the present work we will not address a variety of 
rather important aspects of the electron-phonon inter- 
action in cuprates. In calculating the electron-phonon 
coupling constants, we will completely ignore the corre- 
lated nature of the involved charge carriers, even though 
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this is expected to be essential to deduce reliable quanti- 
tative information about the strength of this interaction 
and its relative importance for different phonon modes. 
Based on calculations by Zeyher and KulidlZI, we can at 
least show that the latter aspect does not oppose to the 
conclusions of our results. Furthermore, lattice anhar- 
monicities will also be neglected, since we believe that 
within the presented framework, where the mutual influ- 
ence of electron- phonon and spin fluctuation interaction 
is investigated within an Eliashberg theory, anharmonic- 
ities, if relevant for the cuprates at all, can hardly affect 
our result. These rather drastic approximations are nec- 
essary to enable us to draw quantitative, material spe- 
cific conclusions about the role of phonons within a spin 
fluctuation scenario. The main complication in this con- 
text is that, in distinction to the electronic system, the 
phononic degrees of freedom behave completely three di- 
mensional, and there is no obvious restriction to a few 
relevant phonon modes. 

II. THEORY 

The coupled electron-phonon system is characterized 
by the Hamiltonian: 

H = ^£k4,/k,,T + gsf ^2 C k,<T^' C k+q.<T' • S_ q 
kcr kqcrcr' 

+ 9A(q) C L C k+q,<T (*>q,A+&q,.\) 

k,qAo" 

qA 

Here, (£ is the electron creation operator of the hy- 
bridized Cu-3d :E 2_j / 2 and 0-2p Xi y orbitals with energy rSt 
describing the fermionic low energy degrees of freedomlla. 
The spin fluctuations are treated in terms of the interac- 
tion of the electron spin with a spin-field Sq characterized 
by a spin susceptibility x(q, a;). Note, alternatively one 
could also use a Hubbard like Hamiltonian and deter- 
mine x(q, us) diagrammatical^, e.g. within the fluctua- 
tion exchange approximatior£3, which would not change 
the conclusions of this paper. The present method has 
the advantage that x(q, lj) can be chosen in agreement 
with the experimentally observed spin susceptibility of 
YBa2Cu3C>7. In Eq. (|l|) o qA is the creation operator of 
a phonon in the vibrational branch A with energy tOq^x 
and gx (q) is the electron-phonon coupling constant of the 
A-th phonon branch with the above Cu-0 hybrid states. 
For the electronic band dispersion we use 

e(k) = — 2t (cos k x + cos k y ) — At' cos k x cos k y , (2) 

with the nearest neighbor hopping t — 250 meV and 
the next nearest neighbor hopping t' = — 0.45f, which 
reproduce the bandwidth and Fermi surface shape of 
YBa 2 Cu 3 7 . 



The short ranged antiferromagnetic spin fluctuations 
are described within the nearly antiferromagnetic fermi 
liquid approach by the dynamical spin susceptibility^: 

Imx(q, cj) = Im ( - XQ — — I . (3) 

VI + ^(Q - qr - iLu/uj s f J 

Here, values for the antiferromagnetic correlation length 
£ = 2.3a, the spin fluctuation energy_a; s f = lAmeV 
and XQ = 75meV are chosen as in Ref.tj, which are de- 
duced from the analysis of NMR and NQR experiments 
of YBa2Cu3C>7. The resulting pairing interaction of the 
spin fluctuation interaction is given by: 

V^( W )= 9s 2 /X (q,w). (4) 

The electron-spin fluctuation coupling constant g s f — 
0.68eV for an electron concentration of n e = 0.8 is cho- 
sen as in Rem, which was shown to yield T c = 90-fT 
and quantitative agreement with transport experiments, 
if one .assumes that these are entirely due to spin fluctu- 
ations^. 

The phonons are described within a shell model, where 
the dynamical matrix is obtained from the derivative of 
the lattice potential, which consists of Coulomb poten- 
tials with core and shell charges, Born-Mayer repulsive 
potentials and core-shell polarizabilities. Details of the 
calculation are the same as in Ref.El For the description 
of YBa2Cii307 we use parameters for the core and shell 
charges, Bor-n.-Mayer parameters and polarizabilities as 
given in Ref.E3. The shell model was shown to give excel- 
lent agreement with experimentally determined phonon 
dispersion curves for several cuprate systems and is based 
on a rather small set of physically well motivated parame- 
ters. The corresponding electron-phonon coupling is cal- 
culated using a model recently proposed by ZeyhertJ, 
which uses rigidly displaced ionic potentials screened by 
a dielectric constant and by holes within the C11O2- 
planes. We have modified this model to allow for the 
explicit consideration of local field corrections taking the 
relative position of the screening CitC^-planes and the 
vibrating ions explicitly into account. The coupling of 
the phonon-mode A to an electron localized at the site k 
within the primitive cell is given by: 

t \ tt^ Z K ,e 2 / h /e are / A (q) ( 2 ) 

-^<euo)) . ( 5) 

M K and Z K are the mass and the charge of the ion with 
basis index k, v c is the volume of the primitive cell, lux(q) 
and e aK \ are the vibration frequency and the normalized 
eigenvector of the phonon mode A with wave vector q 
and a is a Cartesian index. $^K?a denotes the contribu- 
tion of the of the K'-ion to the total Coulomb potential 
at lattice site k, which is caused by the vibration of the 
ion k' and of the ion k, respectively: 
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Gi,G2 are reciprocal lattice vectors and t?(k) = 
denotes the bare Coulomb interaction. The inverse di- 
electric function e _1 will be discussed in the appendix. 

For simplicity we assume, that the Cooper pair form- 
ing electrons are mainly located at the in-plane Cu-sites. 
Therefore only the coupling g K \ where k refers to the 
Cu(II)-site will be considered. This index k can be sup- 
pressed in the following. The electron-phonon induced 
pairing interaction is given by: 
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Here, the summation with respect to q z takes into ac- 
count that the pairing interaction can also be gener- 
ated by phonons propagating perpendicular to the Cu02 
planes, even if the paired electrons are within these 
planes. 

The total electron-phonon coupling strength may then 
be characterized e.g. via the dimensionless parameter 
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where the McMillan function a 2 F is defined by: 
1 
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| ff (k,k',A)| 2 <5(c-c A (k-k')). 
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Having discussed the spin fluctuation and electron- 
phonon induced interactions, we investigate the super- 
conducting state within the strong coupling Eliashberg 
theory, where as usual, the self energy is expanded with 
respect to unitary and Pauli matrices in Nambu space t % : 

E k (^n) = Y k [ioj n )f° + A k (^ n )f 3 + $ k (iw„)f 1 . (10) 

Y k (iu) n ) = iuj n (1 — Z k (iu) n )) and X k (ioj n ) refer to the 
self energy renormalizations which occur also in the nor- 
mal state and $ k (ia; n ) = A k (ia> n )Z k (iw„) to the anoma- 
lous self energy which is nonzero only in the supercon- 
ducting state. uj u = (2n + l)ir/{3 are the fcrmionic Mat- 
subara frequencies, with (3 = 1/T. The self energies are 
determined by the Eliashberg equations: 
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Vk,k'(* w n) = V k [ k ,(ico n ) - V^_ k ,(iu> n ) , 

and the denominator is given by 

D k /(iuv) = {iu n >Z k ,(iuj n ,)) 2 

-(e k / + X k/ (iu} n ,)) 2 - $ k '(^n' 



(12) 
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After analytical continuation to the real frequency axis, 
iuj n = i(2n + l)ir / j3 — > lo + i0 + , this set of coupled equa- 
tions is solved-self consistently using the numerical frame- 
work of Ref.Eil. We do not restrict the momentum sum- 
mation to the Fermi surface but take the entire BZ into 
account. Since both pairing interactions V t '* and V ph are 
chosen to reproduce the experimental spin-susceptibility 
and phonon-dispersion, we do not take additional renor- 
malizations of Vq(iuj n ) and Vq(iu> n ) into account. 



III. RESULTS 



First, we will demonstrate that our treatment of the 
spin fluctuations yields indeed a superconducting state 



with cL 



symmetry and a reasonable critical tem- 



perature. The effective interaction due to spin fluctu- 
ations is repulsive for all momenta and is peaked at 
the antiferromagnetic wave vector Q = (tt/o, 7r/a), as 
can be seen in Fig. [[} In agreement with Refij, the 
self-consistent solution of the Eliashberg equations yields 
d x 2_ y 2 -superconductivity, as can be deduced from the 
momentum dependence of the off-diagonal self-energy 
in Fig. ||(a), which exhibits the typical form: <I> ~ 
(cos k x — cos k y ) . 

Variation of the temperature shows that the super- 
conducting gap vanishes around 11 OK. Although this is 
slightly larger than the experimental critical temperature 
of Y ' Ba^Cu^O^ , we decided to accept this still reason- 
able value for T c , since the spin fluctuation interaction 
has alread#,been chosen in agreement with transport ex- 
perimentsc3 and manipulation of the electron-spin fluc- 
tuation coupling constants would introduce new and ar- 
bitrary parameters. 

In the following we investigate the influence of the 
electron-phonon interaction on the spin-fluctuation in- 
duced superconducting state. Solving the Eliashberg 
equations yields d^.^-superconductivity even in the 
presence of phonons, as can be seen from the momentum 
dependence of the gap-function in Fig. |^. For obtaining 
some information about the isotope effect, we have per- 
formed Eliashberg calculations for two different oxygen- 
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FIG. 1. Spin fluctuation interaction at T=70K: (a) mo- 
mentum dependence of the spin fluctuation interaction for 
frequency to = 0, (b) frequency dependence of the on-site in- 
teraction Vbo and the nearest neighbor interaction Vio- 

isotopes 0(16) and 0(18). For this the whole electron- 
phonon interaction, i.e. the phonon dispersions and the 
electron-phonon coupling constants, has been calculated 
separately for the two oxygen masses. In Fig. || the mo- 
mentum averaged density of states at a temperature of 
T = 80A' is shown for the two oxygen isotopes 0(16) 
and 0(18). For comparison the density of states for only 
spin fluctuation induced superconductivity is shown as 
well. Two important things can be learned from Fig. 0: 
the superconducting gap is reduced dramatically in the 
presence of phonons, but nevertheless the isotope effect 
is positive. This can be infered from the fact that there 
still remains a small superconducting gap in the density 
of states of the oxygen isotope 0(16) at the temperature 
chosen in Fig. || (T=80K), whereas for 0(18) the super- 
conducting gap has already completely disappeared. 

Note that similar results -have been obtained by H.- 
B. Schiittler and C.-H. PaoEi They also hnd that the 
electron-phonon interaction reduces the critical temper- 
ature of spin fluctuation-induced superconductivity con- 
siderably, whereas the isotope exponent remains small 
but negative. However in their calculation this is due 
to the local nature of the electron-phonon interaction 
whereas the focus of our work is to investigate the role 
of the k-dependence of this interaction. 

In order to obtain an estimate for the critical tem- 



perature of the two oxygen isotopes, we have performed 
succesive calculations in the temperature range between 
78K and 83K. Since the self-consistent solution of the 
Eliashberg equations gives rise to a finite single parti- 
cle scattering rate in the gap once the temperature is 
nonzero we choose the disappearance of the dip struc- 
ture in the DOS as criterion for T c . In this way we ob- 
tain for the oxygen isotope 0(16) a critical temperature 
of T c = 82K and for the isotope 0(18) T c = 80K. Using 
these values for the calculaion of the isotope exponent 
a = — ~ — t^s^j we obtain an approximate 

value of a sa 0.2, which is only slightly larger than exper- 
imentally found in optimally doped cuprates. 

For a better understanding of the phononic contribu- 
tion to the pairing interaction and the reason for the 
positive isotope effect, the phonons will be analyzed sepa- 
rately now. The phonon density-of-states and the McMil- 
lan function a 2 F are shown in Fig. ||. For total electron- 
phonon coupling strength we obtain A = 0.47, which is. 
in good agreement with estimates of Xtot = 0.4— 0.6E§E£l 
from experimental determination of frequency shifts and 
line widths. Nevertheless, the calculated A = 0.47 should 
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FIG. 2. Off diagonal self-energy $fc(a») for (a) only 
spin fluctuation induced superconductivity and (b) both 
phonon and spin fluctuation induced superconductivi- 
ty. In both cases the superconducting state is char- 
acterized by d x 2_ y 2 -symmetry order parameter, since 
<&k ~ (cosk x — cosky). Note, $fc(u) is reduced due to the 
electron-phonon interaction. 
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we show the on site interaction Vr 
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FIG. 3. Momentum avaraged electronic density of states 
at T=80K for a superconducting state induced only by spin 
fluctuations and by spin fluctuation and electron-phonon in- 
teraction together for two different different oxygen isotopes 
0(16) and 0(18). 

lead to a considerable contribution to the resistivity, 
which within the spin fluctuation model is believed to 
be dominated by purely electronic scattering. This indi- 
cates that role and relative weight of these two scattering 
channels for transport phenomena is not yet completely 
understood. We checked that modifications of A due to 
backscattering vertex-corrections do not reduce A consid- 
erably. 

In Fig. ||(a), we show the momentum dependence of 
the zero frequency effective phonon mediated electron- 
electron interaction of Eq. (0) . The interaction is attrac- 
tive for all momenta and most important for the present 
study peaked for small momentum transfer. In Fig. |5](b) 
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FIG. 4. Phonon density of states for Y BaiCuiOi ob- 
tained from a shell-model calculation using parameters as in 
Ref. [22], and the McMillan function a 2 F(co). The calcula- 
tion of the DOS and a 2 F(ui) were performed with a channel 
width of 0.2 meV. 
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l £ V(k X ,ky) , 



(15) 



which represents the pairing interaction in the case 
of isotropic s-wave superconductivity, and the nearest 
neighbor interaction Viq: 



'io 



— ^ V(k x ,k y )cosk x , 



(16) 



which represents the pairing interaction in the case of 
c?a;2_ a 2-symmetry superconductivitycJ. Both Voo and Vio 
are attractive in some frequency range, but the on site 
interaction is much stronger and the corresponding at- 
tractive frequency range is slightly larger. Therefore, 
phonons alone will give rise to s-wave superconductiv- 
ity. This is natural for a phononic mechanism, since 
the electron-phonon interaction is attractive for all mo- 
menta and thus the momentum averaged interaction Voo 
is dominant. The self-consistent solution of the Eliash- 
berg equations confirms this. Taking only the electron- 
phonon interaction into account, we obtain superconduc- 
tivity with anisotropic s-wave symmetry without nodes 




Frequency [meV] 

FIG. 5. Electron-phonon interaction at T=70K and 
A p h = 0.47: (a) momentum dependence of the elec- 
tron-phonon interaction Vj. for frequency uj = 0, (b) frequency 
dependence of the on-site pairing interaction Voo and the near- 
est neighbour pairing interaction Vio. 
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T=45K A. nh =0.96 
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FIG. 6. Off diagonal self-energy for only phonon induced 
superconductivity. The momentum dependence exhibits 
clearly an anisotropic s-wave symmetry. 

as can be seen from the momentum dependence of the off- 
diagonal self energy in Fig. [| For obtaining this result at 
the numerically accessible temperatures T > 30 K we had 
to double the total electron-phonon coupling strength ar- 
tificially. But we expect that enhancing the coupling 
strength should solely increase the critical temperature 
and not change the symmetry of the superconducting 
state. Even for this artificially large coupling strength 
of A = 0.96 superconductivity already vanishes between 
about T = 45 — 50K. Hence, for realistic coupling 
strengths the critical temperature of phonon induced su- 
perconductivity is much smaller than the critical tem- 
perature of spin fluctuation induced superconductivity. 
Note that in this case the repulsive Coulomb interaction 
was not taken into account, which will suppress the criti- 
cal temperature of the s-wave superconducting state even 
more. 

Most interesting, however, is the result that all phonon 
modes together can generate an attractive interaction be- 
tween nearest neighbor sites. Thus the consideration of 
electron-phonon coupling increases the pairing interac- 
tion in the d x 2_ y 2-wave channel, which is a necessary 
prerequisit for a positive isotope effect. The reason for 
the resulting decrease of the superconducting gap due 
to the presence of phonons is the large on-site interac- 
tion Vqo, which enhances the scattering rate dramati- 
cally. Therefore, concerning the sign of the isotope ef- 
fect, it is plausible, that the isotope effect is positive 
as long as the positive contribution from the increased 
pairing interaction Vio (where the larger cut-off corre- 
sponds to the smaller oxygen mass) outweighs the nega- 
tive contribution from the increased scattering rate due 
to larger on-site interaction Voo- Also it is clear that 
including both spin fluctuation and electron-phonon in- 
teraction still yields d x 2_ y 2 -superconductivity, since the 
consideration of electron-phonon coupling even increases 
the d x 2_y 2 -pairing interaction, whereas the spin fluctua- 
tion contribution to VoO) the s-wave pairing interaction, is 
strongly repulsive, as can be seen from Fig. [jj and there- 
fore reduces the s-wave pairing interaction dramatically. 



We will focus now on the question, why the near- 
est neighbor electron-phonon interaction is attractive at 
all. The attractiveness of Vio can be traced back to the 
fact that the phononic interaction is peaked for small 
momentum transfer, since multiplication with cos k x in 
Eq. (|ll]) gives more weight to small momenta. Ob- 
viously, besides the shape of the Fermi surface, which 
is taken from ARPES experiments, the momentum de- 
pendence of the electron-phonon interaction is the key 
feature, which controls possible contributions to differ- 
ent pairing symmetries. The momentum dependence of 
the electron-phonon interaction consists of two contribu- 
tions, as can be seen by inspection of Eq. (^); namely the 
momentum dependence of the electron-phonon coupling 
constants and the momentum dependence of the phonon- 
Green's function. Because the momentum dependence 
of the phonon-Green's function is directly related to the 
dispersion of the phonon spectrum, which is relatively 
weak in the case of dominating optical phonon modes, 
we find that the momentum dependence of the electron- 
phonon coupling constants is the decisive factor. There- 
fore, phonon modes with coupling strength largest for 
small momentum transfer can enhance the d x 2_ y 2 pair- 
ing interaction. However, phonon modes with coupling 
strongest for large momentum transfer reduce the d x 2_ y 2 
pairing interaction. 

In order to find out which phonon modes support the 
d 3 .2„j / 2-pairing interaction, the vibrations of single atoms 
(per unit cell) will be investigated. For illustration of the 
structure of YBa2Cu30y we refer to Fig. fj]. The great- 
est coupling strength is caused by the vibration of the 
in-plane oxygens O(II), with vibration amplitudes along 
the Cu-0 bonding axis and along the z-axis, i.e. breath- 
ing mode like and buckling mode like vibrations, and by 
the vibrations of the apex oxygen 0(1) in z-direction. For 
illustration of the displacement patterns of the breathing 
mode and the buckling mode we refer to Fig. ||. In Fig. ^ 
the momentum dependence along the T — ► M -direction 
of the corresponding electron-phonon coupling is shown. 




FIG. 7. Structure of YBa 2 Cu 3 7 . 
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The in-plane oxygen 0(11, III) vibrations are quite well 
described by the momentum dependence g(k) ~ (sin k x + 
sinfe„) and g(k) ~ (cosfca, + cosk y ), which is the re- 
sult for the breathing and buckling mode, respectively, if 
only nearest neighbor interaction is considered. Thus the 
breathing mode reduces d x i -superconductivity, be- 
cause the momentum dependence of its coupling strength 
is peaked for large momentum transfer. On the other 
hand the coupling strength of the buckling mode is 
peaked for small momentum transfer and therefore the 
buckling mode supports (4,2 _ a 2 -wave superconductivity 
in so far as the d x 2_ y 2 -pairing interaction is increased by 
the buckling mode. The small wiggles in Fig. |^ are caused 
by wiggles in the momentum dependence of the free sus- 
ceptibility which determines the screening proper- 
ties of the Cu02-planes (see appendix). In the case of 
the apex oxygen vibrations, the momentum dependence 
of the coupling is different for the two apex ions in the 
primitive cell: the closer apex oxygen causes a more or 
less momentum independent coupling, whereas the cou- 
pling of the further apex oxygen is stronger for small 
momentum transfer similar to the buckling mode like 
vibrations. In respect to vibrations of the other ions, 
which are not shown in Fig. we find that the coupling 
of vibrations in z-direction is generally peaked for small 
momentum transfer. The vibration of the oxygen ions 
is the most dominant contribution to the total coupling 
strength, whereas the contribution of yttrium and barium 
is very small due to screening in the Cu02-planes. Alto- 
gether all vibrations in z-direction, especially the buck- 
ling modes, increase the nearest neighbor interaction Vio- 
Only the breathing mode and the vibration of the chain 
oxygen O(IV) along the CwO-bonding axis reduce Vio- 
The other modes are of no relevance with respect to Vio- 

Note that within our approach the coupling strength 
of the buckling mode is about twice as large as the cou- 
pling strength of the breathing mode. The situation is 
opposite, if the same calculations are performed with the 
original Zeyher-model. The reason for this is to be found 
in the different treatment of screening. In Zeyher's model 
the screening planes are smeared over the whole elemen- 
tary cell, whereas in our approach the exact position of 
the Cu02-planes with respect to the vibrating ion and 
the Cooper-pairing hole is considered. It is evident that 
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r M 
FIG. 8. Displacement patterns of the buckling mode and 
of the breathing mode at the F-point and the M-point of the 
BZ. 
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FIG. 9. Electron-phonon coupling constant for the vibra- 
tion of the in-plane oxygen O(II) (a) in x-direction, (b) in 
z-direction and (c) for the vibration of the apex-oxygen 0(1) 
in z-direction. Since every hole is located just in one of the 
planes, the coupling is different for the vibrations of the oxy- 
gens of the two planes. The coupling of the oxygen in the 
same plane is denoted by a solid line and the coupling of the 
oxygen of the other plane by a dashed line. Note, in the case 
of the apex-oxygen, this refers to the closer and the further 
apex-oxygen. 

the screening of the breathing mode is better, if one ac- 
counts for the fact that both the vibrating oxygen ion 
and the Cooper pair forming hole are located within the 
same plane. On the other hand the opposite is true for 
the buckling mode. Because the vibration of the oxygen- 
ion in z-direction is perpendicular to the Cit02-plane, it 
cannot be screened within the same plane, but only by 
charges induced in more distant planes. Therefore, the 
buckling mode is less screened, if one takes the relative 
location of the vibrating oxygen-ion, the Cooper-pairing 
hole, and the screening plane explicitly into account. Ad- 
ditionally, it has to be pointed out that in our approach 
(as in the Zeyher-model) only the electrostatic part of the 
electron-phonon coupling is considered, whereas all con- 
tributions arising from changes in the overlap of neigh- 
boring electronic wave functions have been neglected. It 
is to be expected that these neglected effects are strongest 
for the breathing mode. This means that the coupling 
strength of the breathing mode could be somewhat un- 
derestimated in our approach. In this respect it is inter- 
esting to note that recently the electron-phonon coupling 
was calculated by O. K. Andersen et alEOl and O. Jepsen 
et alEB within a tight-binding model derived from LDA 
linear-response calculations. They find a total electron- 
phonon coupling strength of A s f=a 0.4, which is similar to 
our findings, and a medium d x 2 _ y i -wave-coupling strenth 
of Xd ~ 0.3, which is even much larger than our value. 
Most interesting, they also find a very dominant cou- 
pling of the buckling mode. Although their approach is 
complementary to ours, their results are quite similar, 
which supports our results concerning the relative cou- 
pling strengths of the breathing and the buckling mode. 
Anyway we expect that a more dominant breathing mode 
would reduce Vio somewhat, leaving the other results 
qualitatively unchanged. 

Our numerical results, obtained within a strong cou- 
pling Eliashberg theory, show that the electron-phonon 
interaction reduces the critical temperature, even though 
it increases the d-wave pairing interaction and the isotope 
effect is positive. In order to investigate these aspects in 
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a qualitative, but more transparent way, we use a sim- 
ple weak coupling model, where the pairing interactions 
caused by spin fluctuations and phonons are approxi- 
mated by frequency independent coupling constants A^ 
and Ap h with different cut-offs w s f and ui p h (Fig. |l0|). In- 
troduction of the renormalization constants Z s f = l+X s j 
and Z s f—ph = 1 + A s / + X p h accounts for the reduction of 
the spectral weight of the Cooper pair forming holes by 
scattering effects due to spin fluctuations and phonons. 
For the renormalization factors the same cut-off is used as 
for the pairing interaction. The coupling constants A are 
defined in Eq. ||. The coupling constant X d characterize 
the d x 2_ y 2 -pairing interaction: 



A" 



E kk ^k)^kQiw(k,k',o)<^ 

Ekk'<^k)£(ek')^ 



(17) 



with ip^ — cos k x — cos k y . These approximations result in 
the following form of the gap-function A k = Ao?/>^ with 
A = Ai if |e k | < ui ph and A = A 2 if uj ph < |e k | < u sf , 
respectively. A 1; A 2 are obtained by solving the gap 



equations: 



A, 



A ?/ + ^7V-2/(o Al 



-"sf-ph 







* f de l^IM A2 



and 



A 2 = 



V 



Sf-ph 



C2 



ei 



£2 



(18) 



(19) 




FIG. 10. Modelling of the spin flucutuation interaction 
and the electron-phonon interaction by frequency indepen- 



dent coupling constants. Here, V a f = Afy and V p h 



denote the d a .2_ H 2-type contribution of the corresponding in- 
teractions, respectively with cut off u a f and lu p h ■ 
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FIG. 11. Results for A and T c using the simplified inter- 
action model taking u 3 f = 200meV, Xf* = 1.1, X s f = 3.2, 
\ d ph = 0.3, \ ph = 0.6. 

with 



z 2 
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A? 
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(20) 



In the following we will study the influence of the size 
of the phonon cut-off us p h on the superconducting prop- 
erties by means of this simple model. In Fig. |ll](a) the 
gap-functions Ai and A 2 corresponding to the regions I 
(|e k | < Lu ph ) and II (uj ph < |e k | < u s f) at T — OK are 
plotted and in |ll|(b) the corresponding critical tempera- 
ture is shown. The coupling constants A s /, Xjp X P h, X& h 
are chosen in agreement with the previous results. Only 
Xp h is chosen slightly larger for illustrational purposes. 
Since the spin fluctuation interaction is decreasing slowly 
over a large frequency range, as can be seen by inspec- 
tion of Fig. H[ it is hard to estimate the spin fluctuation 
cut-off uj s f. In this calculation uj s f = 20077161^ was cho- 
sen as a cut-off. Two things can be observed in Fig. |Tl] : 
First, the gap-function Ai (region I) is generally larger 
than the gap-function A 2 (region II) due an increased 
pairing interaction caused by phonons. Furthermore the 
gap-functions Ai and A 2 and the critical temperature 
T c show a peculiar cjp/j-dependence: for small w p h they 
are decreasing as a function of u) p h and after a certain 
critical cut-off frequency they are increasing again. Ex- 
amination of gap-equation Eq. (|l^) shows that there are 
two contributions, which determine the size of the gap- 
function Ai: the interaction of charge carriers(I) among 
themselves and the interaction of charge carriers(I) with 
charge carriers(II). In an analogous manner there are 
two contributions, which determine the size of the gap- 
function A 2 . Only one of these (four) contributions, i.e 
both electrons in region I, can profit from the phononic 
part of the pairing interaction. On the other hand two of 
these (four) contributions, each of which result from the 
interaction of the respective electrons with electrons of 
region I, are reduced, because the quasi particle spectral 
weight in region I is diminished by the additional scat- 
tering on phonons. Altogether only the gap-function Ai 
can profit from the phononic contribution to the pairing 
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FIG. 12. Dependence of the isotope exponent a on the 
spin fluctuation d c 2_ !/ 2-coupling strength and the total 
spin fluctuation coupling strength, respectively. These re- 
sults were obtained with the simple weak coupling model 
with u a f — 200meV, u) p h = 70meV and additionally (a) 
A s / = 2.5 and (b) Af/ = 1.3, respectively. The total elec- 
tron-phonon coupling strength was taken as X p h = 1, but 
the rf a .2_ J/ 2-electron-phonon coupling strength was varied as 
indicated in the figure. 

interaction, whereas both gap-functions Ai and A2 no- 
tice the destructive effect due to scattering on phonons. 
This is the reason why for small uj p h both of the gap- 
functions Ai and A2 are decreasing as a function of oj p h- 
The larger region I the more electrons are scattered by 
phonons. Only when region I has become large enough 
and uiph exceeds a critical cut-off frequency, the positive 
effect of increased pairing interaction between Cooper 
pairs in region I outweighs the negative effect of addi- 
tional scattering and thus both gap-functions Ai and A2 
start increasing as a function of u) p % for larger frequencies. 

The Wp/j-dependence of T c shown in Fig. pT|(b) displays 
an analogous behavior. For very small Lo p h the critical 
temperature is decreasing, and after passing through a 
minimum, T c is increasing again yielding a positive iso- 
tope effect. The critical temperature itself will, for real- 
istic phonon frequencies uj p h, be smaller than the critical 
temperature of only spin fluctuation induced supercon- 
ductivity. 

From Fig. |ll](b) we can directly extract the isotope 
exponent. For very small phonon frequencies the criti- 
cal temperature is decreasing for increasing u p h- Hence, 
the isotope exponent a is negative in this region. For 
larger uj p h, after passing through a minimum, T c is in- 
creasing with u)ph- This causes a positive isotope ex- 
ponent. The magnitude of a depends on the relative 
size of the coupling strengths of the electron-phonon and 
spin fluctuation interaction and the corresponding scat- 

V 



tering rates. In Fig. [12j the isotope exponent vs Af f and 



A s / is shown for different values of \^ h . For decreasing 
spin fluctuation contribution Af* to the pairing interac- 
tion the isotope exponent is increasing strongly. This is 
obvious, because the electron-phonon coupling becomes 
the origin for superconductivity, a situation which seems 
unlikely for cuprates due to the large phononic coupling 
constants necessary in this case. 

The isotope exponent of the present work is al- 



ways small in magnitude and positive. For the system 
YBa2Cu307, under consideration, this is in qualitative 
agreement with the experimental observation. Whether 
the increase of a for non-optimal doping concentration 
is solely due to the decreasing strength of the leading 
electronic pairing interaction or of more subtile nature is 
beyond the scope of the present paper. 



IV. SUMMARY 



In summary, we find that spin fluctuation induced 
d E 2_ y 2-superconductivity is robust against the existence 
of electron-phonon coupling. The superconducting tran- 
sition temperature is reduced by the interaction with 
phonons due to an enhancement of the electron scat- 
tering rate. Furthermore, we find that the momentum 
dependence of the electron-phonon interaction increases 
the d x 2 _ y 2 -pairing interaction slightly. This is important 
for getting a positive isotope effect. Our results indicate 
that the isotope effect is small for optimal doping. If 
one assumes that the spin-fluctuation coupling strength 
varies with the doping concentration, this would imply a 
strong doping dependence of the isotope effect. 

A major approximation of the present work is the ne- 
glect of the effect of strong electronic correlation on the 
electron-phonon coupling constants. As shown by Licht- 
enstein and KulicH strong electronic correlation cause 
a suppression of charge fluctuations at small distances. 
Therefore, we expect g\(k) to be reduced for zone bound- 
ary phonons. This would decrease X ph and thus the large 
scattering rate due to phonons. For the same reason, we 
expect X ph to be slightly enhanced by electronic correla- 
tions. 

Our results demonstrate that the existance and sign of 
the isotope effect can be understood within a spin fluc- 
tuation induced d-wave state. Due to the d-wave nature 
of the superconducting state, it is particularly impor- 
tant to obtain a proper description of the momentum 
dependence of the electron-phonon interaction. It was 
shown that out-of-plane vibrations are essential for the 
d-wave state and that one cannot focus solely on one 
particular phononic mode. Finally, the strong electronic 
correlations inhibit the occurrence of a phonon induced 
superconducting state with s-wave symmetry. A phonon 
induced d-wave state can be excluded because of the in- 
sufficient strength of this subdominant pairing interac- 
tion. 
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APPENDIX A: DERIVATION OF THE INVERSE 
DIELECTRIC FUNCTION 

In this appendix the treatment of the local field cor- 
rections arising from the position of the screening Cu02- 
planes relatively to the vibrating ion will be discussed. 
If spatial inhomogenity is taken into account, the dielec- 
tric function becomes a matrix in the reciprocal lattice 
vectors G, G': 

e(k + G,k + G') = eoo(5GG<- 

w(k + G)x(k + G,k + G') . (Al) 

Here, x is the irreducible polarization of the density- 
density correlation function x, which characterizes the 
screening properties of the system. The inverse dielectric 
function is defined by: 

e(k + G, k + G")e _1 (k + G", k + G') = <5 GG < . 



G" 



(A2) 



This can be expressed as: 



e-^k + G.k + G') = — (<S, 



GG' 



-v(k + G)x(k + G,k + G')) 



(A3) 



The susceptibility X has to be determined by solving the 
following integral equation: 

eoo x(k + G, k + G') = x(k + G, k + G') + 

X(k + G, k + G>(k + G")x(k + G", k + G'). 



G" 



(A4) 



In the cuprates the most dominant screening is caused by 
mobile charges within the Cu02-planes. In YBa2Cu307 
there are two planes per unit cell. Therefore, the sus- 
ceptibility x is a 2 x 2-matrix, describing density-density 
correlations within the same plane and between both of 
the planes. The polarization x on the other hand reduces 
to a diagonal matrix, because the particles are assumed 
to be confined to the planes. In real space the polariza- 
tion is given by: 

X(i"i, r 2 ) = X 2D (xi - x 2 )5(z 1 - z 2 ) ^ 8{z l - z jP ) , 

(A5) 

where Zj@ = jc + dp are the z-coordinates of the CuOi- 
planes, with j being the index for the primitive cell, j3 the 
index for the Cu(92-plane within the primitive cell and c 
the height of the primitive cell. Fourier transformation 
yields: 



X (k + G, q + Q, q + Q') = -x 2D (^ + G 

c 



Here, k, G refer to the in-plane wave vectors, whereas 
q, Q represent the wave vectors in z-direction. In RPA 
X 2D reduces to the noninteracting particle-hole bubble 
of the holes in the CitCVplanes. x is diagonal with re- 
spect to the in-plane wave vectors and the Q-dependencc 
is simple enough. This permits an anal ytic treatment of 
the Q-summations. Therefore, Eq. (A4) can be inverted 
easily and the inverse dielectric function can be obtained 
in this way. There remains only the question of how to 
choose the z-coordinate dp of the screening planes, which 
are not unambiguous due to the internal buckling of the 
Cu02-planes. In our calculations we assumed that the 
screening plane should follow the buckling of the CuOi- 
plane. 
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